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We calculate a novel "magnetic contribution" to the dilepton spectrum in heavy-ion colli- 
sions arising from interaction of relativistic quarks with intense magnetic field. Synchrotron 
radiation by quarks, which can be approximated by the equivalent photon flux, is followed 
by dilepton decay of photons in intense magnetic field. We argue that "magnetic contribu- 
tion" dominates the dilepton spectrum at low lepton energies, whereas a conventional photon 
dilepton decay dominates at higher lepton energies. 

1. INTRODUCTION 

Electromagnetic radiation that accompanies any relativistic heavy ion collision weakly interacts 
with hot nuclear matter. However, it strongly interacts with highly intense magnetic field that is 
frozen into the nuclear matter [1-4]. Therefore it bears witness to the magnetic field existence and 
provides a rare opportunity to experimentally study its properties. In practice, electromagnetic 
radiation caused by magnetic field is masked by electromagnetic radiation of quark-gluon plasma 
(QGP). It is not an easy problem to disentangle the two contributions. A good theoretical control 
over both is required to achieve this goal. 

In the present article wc address the problem of dilepton production in heavy-ion collisions due 
to magnetic field. Our work is partly motivated by the resent experimental results on dileptons and 
photons produced in heavy-ion collisions, which challenge the prevailing theories of electromagnetic 
processes in high energy nuclear physics that neglect strong electromagnetic interactions. It was 
observed in [5, 6] that there is a significant underestimation of electromagnetic spectra by theoretical 
models in the low momentum region. This indicates that there are additional contributions that 
have not been taken into account. We dub the additional lepton production due to magnetic field 
as the "magnetic contribution" as opposed to the "conventional contributions" . 

Dilepton production in magnetic field proceeds in two stages. Firstly, quark (or antiquark) 
is produced in a heavy-ion collision. We will refer to the quark (or antiquark) distribution at 
this stage as the initial quark distribution. Secondly, it radiates the lepton-anti-lepton pair in a 
processes mediated by a virtual photon as depicted in Fig. 1 and Fig. 2. Double lines indicate quark 
and lepton propagators in magnetic field. In Fig. 1 both quark and lepton move in magnetic field. 
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FIG. 1: Dilepton production by quark in external magnetic field, mediated by a virtual photon. Double 
lines denote fermion propagator in magnetic field. Components of four-vectors: p = {£,p}, k = {ui, fe}, and 
P± = {£±,P±}- 

whereas in Fig. 2 only one of them does. The final result is sum over all these processes. However, 
in this article we will focus on the square of the amplitude shown in Fig. 1, because we anticipate 
that it is most sensitive to magnetic field. Other contributions will be considered elsewhere. 



The initial quark distribution is made up of soft and hard quarks and antiquarks. Soft quarks 
have typical energies e ~ T, where T is the plasma temperature, and are a part of QGP. The 
exact mechanism of their production and equilibration is not fully understood and is actually not 
essential for our arguments. What is important is an observation that according to the state-of- 
the-art phenomenology, the QGP equilibration happens over very short time on the order of 1/Qs, 
where Qs is the saturation momentum, see e.g. [7]. Hard quarks in the central rapidity region 
y = have typical energies e ^ Qs ^ T and their production mechanism has been thoroughly 
investigated [8-14]. In this article we focus on lepton energies e+ in the interval T < e-^- < Qs- In 
this region, spectrum of soft quarks falls off exponentially, whereas the spectrum of hard quarks is 
only logarithmic. Moreover, the number of soft quarks in plasma of volume V is on the order of 
T^V. At early times V ^ S /Qs where S is the cross-sectional area of the ion overlap region. On 
the other hand, the number of hard quarks is on the order of SQ^s (see Sec. 4), which is much larger 
than the number of soft quarks. Therefore, the contribution of hard quarks to dilepton production 
is dominant. It is this contribution that we discuss in this article. 

The paper is organized as follows. In Sec. 2 we employ the Weiszacker- Williams method to 
write the dilepton production rate by a hard quark as a convolution of the real photon decay 
rate with the flux of equivalent photons emitted by a fast quark. In Sec. 3 we compute the rate of 





FIG. 2: Magnetic contributions that are not considered in this article. 



3 



photon dissociation into dilepton as a function of lepton transverse momentum and rapidity. Initial 
quark distribution produced in heavy-ion collisions is computed in Sec. 4 using the quasi-classical 
approximation. In Sec. 5 we apply the developed formalism to calculate the magnetic contribution 
to e+e" pair production at midrapidity at RHIC. We compare dilepton production in magnetic 
field and in vacuum in Sec. 6 and argue that magnetic field contribution dominates at photon 
energies uj < 200 MeV. Finally, we summarize in Sec. 7. 



2. EQUIVALENT PHOTON APPROXIMATION 

Calculation of dilepton spectrum significantly simplifies because light quarks are ultra- 
relativistic £ S> in the center-of-mass frame. This allows us to employ the Weiszacker- Williams 
method to calculate the dilepton production with logarithmic accuracy. According to this method, 
we can relate the cross section of g ^ qi'^i" process to the cross section of photo-production 
7 — >■ t^l~ . The logarithmically enhanced contribution arises from the kinematic region where 
photon virtuality has negligible effect both on photon emission and on dilepton photo-production. 
The dilepton production rate can be written as (see notations in Fig. 1) 



I niUj) ——r- ; dOJ , (1) 



dt dr2_|_d£_^ 

where n{oj) is the flux of equivalent real photons replacing the virtual photon. It is given by [15] 

2z'^a 1 £ 
n(uj) = In / , , (2) 

where Zg is the quark's charge and 



X = ^ JHF^' = ^\pxB\ = ^p,B (3) 



is a boost-invariant parameter. Here pt is quark's momentum component transverse to B (not to 
be confused with p±, which is transverse to the heavy-ion collision axis). Eq. (1) takes into account 
that the spectrum of photons is enhanced at a; ^ £. Angular distribution of photons strongly peaks 
at small angles ~ rrig/e'^. Therefore, the equivalent photon momentum k is approximately coUinear 
with the quark's momentum p. 

It is convenient to introduce dimensionless parameters 6 = a^^'^B/rrig, which is the value of 
magnetic field in units of the Schwinger's critical field, and 7 = e/iriq. In terms of these parameters 
X = 67 sin 0, where 6 is the angle between p and B. In heavy-ion collisions 7 S> 1 and b > 1 
implying that % S> 1. This estimate breaks down only at very small angles 9 < (67)""^ where 
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X becomes small and the photon decay rate is exponentially suppressed [16] and hence can be 
neglected. Integration in (1) runes over frequencies 2m^ <oj< ojm with a;^ satisfying the equation 



1 + (xe/w^ 



a/3 



= e . 



Since uj <^ e and % S> 1 it follows that ex/w » 1. Using this in (4) we get 



(4) 



(5) 



To determine the region of applicability of the equivalent photon approximation, note that the 
logarithmic contribution to the equivalent photon flux comes about only if Um ^ m^. In view of 
(5) we have 



7 mj 
- ^ — - 
b rrii 



(6) 



This condition is satisfied for not too strong fields and for light leptons. In particular, we will apply 
the equivalent photon approximation to calculate the magnetic contribution to electron-positron 
pair production at RHIC in Sec. 5. 



3. PHOTON DISSOCIATION RATE 



Now wc turn to the photon dissociation rate which is convenient to perform in a frame where 
photon's momentum is perpendicular to the magnetic field. We denote such frame as K' and all 
quantities in it bear the prime. The corresponding dilepton rate is given in Sec. 3 A. In Sec. 3 B the 
rate is transformed to frame K where photon moves at an arbitrary angle with respect to magnetic 
field. 



A. JT'-frame: k' ■ B' = 0. 



In reference frame K' photon moves in a plane perpendicular to the magnetic field: k' ■ B' = 0. 
The rate of dilepton photo-production in K' frame reads [16] 



dN, 



a 



dt'd^p' (27r)2w 



dr 



e'e' 



./2 , /2n 2 2 



xexp{z^a;V(l-fe'.< + '^)}, 



(7) 



where e', + £'_ = uJ and the synchrotron frequency is 



^^0 = — 



(8) 
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Integration over the time parameter r can be done explicitly using the following formulas 

b 



/_^cos[6r + ar3]dr = ^Ai(- 



,(3a)V3;' 

cos[6r + aT^]dT = - — / cos[6r + ar^]dr , (10) 

-oo "J^ J —oo 

where Ai(x) is the Airy function. Using (9), (10) into (7) we obtain 



_ -I- 9 + ~ h _ k' . W 

+ \ 



47r 



e' ^ 



xAi^2(l-*'..;)^ji^j j. (11) 

Integration over the lepton direction, given by velocity v^^, is convenient to do in (7) before taking 
the r-integral. It can be performed with the required accuracy along the steps outlined in [16]. 
The main contribution arises from radiation of collinear lepton, i.e. k' ■ v', 1. The result is 



dK 



dt'de'^ a;'2 



ki{^)di + (^^ - xx'V2^ Ai'(x')} , (12) 



where Ai'(a;) is the derivative of the Airy function and 

3 / \ 2/3 



-=^. (13) 



Parameter x is boost-invariant. In an arbitrary frame it reads as follows 



X = -^J-{F^,k-Y = ^\kxB\. (14) 



dN^^g+g- _ aeB' 2{x'^l^ + 1/ x) {x')_^^, ^^^^ 



The total rate is 



dt' mtX 7(4/^)2/3 x'lV4(2;'3/2 _ 4/^)1/2 

B. if-frame: k ■ B ^ 0. 

Now wc need to transform equations (1),(2),(7) to an arbitrary frame K. This is done by 
making a boost in the magnetic field direction. Let z be the collision axis and y be the magnetic 
field direction. We will use the following notations: a and (3 are polar and azimuthal angles of 
photon and quark with respect to z axis (collision axis); 6 and (f) are polar and azimuthal angles 
of photon and quark with respect to y axis (magnetic field direction) . The same symbols with the 
"+" subscript refer to lepton For example, quark momentum reads 

p = £{x sin a cos /3 + y sin a sin /5 + £ cos a) (16) 
= e{x sin cos <p + y cos ^ + £ sin ^ sin <p) , (17) 
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and similarly for other vectors. 

Suppose that K moves with velocity V = Vy with respect to K'. Boost along the magnetic 
field direction does not change the field, i.e B' = B. The Lorentz transformation formulas for the 
quark momentum read 

p'x=Px, Py = lv{Py + Ve) , p'^ =pz, e' = 7y(£ + Vpy) (18) 

Since momenta of photon and quark are approximately coUinear, in K'-frame p'y ^ k'y = 0. This 
implies using (18), (17) that 

y = -cos^ = -sinasin/?, = (1 - y2N-i/2 ^ 1 _ /^qn 

smfc* 

Transformation of quark and photon energies is given by 

e' = £ sing* = — , = — , (20) 

Analogously to (18) we can write the transformation of the positively charged lepton's l'^ momen- 
tum 

p'+x=P+x, p'+y = lv{p+y + V£+) , p\^=p+^, e'j^=-iv{£+ + yp+y) ■ (21) 
Thus, lepton's energy £+ transforms as 



5^+ = 7y£+(l + V cos6+) = 7y£+(l + FsinQ;_|_sin/3+) . 



(22) 



Consider a relativistic invariant 

k-p+ = cje+{l -k-v+) = uj'e'+il - k' ■ v'+) . (23) 

Using (20) and (22) we find 

- , , 1 — fc • u+ 1 — cos 9 cos 0+ — sin sin 0+ cos(^ — 

" 1 + FCOS0+ " l-cos0cos^+ ' ' 



1 — cos a cos — sin a sin a+ cos(/3 — /3+) 
1 — sin Q sin /3 sin aj^ sin /?_|_ 

Finally, employing transformation of the time interval and the solid angle 



(25) 



dt' = —dt , dn' = ^- 1- ^dn+ (26) 



we obtain that in K' frame 



diV _ 1 dN 

dtdn+d£+ ~ 7^(1 + y cos6'+) dt'dQ'^de'^ 



(27) 
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where e', e'j^ and oj' on the right- hand-side should be replaced with the corresponding expressions 
in K frame using (19), (20), (22), (25). 

Experimental data on dilepton production is usually represented in terms of rapidity y+ and 
transverse momentum p-\-±_ - which are convenient if there is no magnetic field - in place of energy 
£+ and polar angle a+. In the ultra-relativistic limit 



P+L = sjvXx +P\y = sin a:+ , y+ = - Intan ^ . (28) 

Inverting equations (28) yields 

£+ = p+_L cosh y+ , sin a+ = — ^ . (29) 

cosh y+ 

Using (29) in (19), (20), (22), (25) we can transform parameters appearing in (11) to A'-frame as 
follows 

sin/3+ sin/3 coshy 

cos6'+ = — , V = — , 7V = , ^ , (30) 

coshy+ coshy ^cosh^ y - sin^ p 

£+ = I (cosh y cosh y+ - sin /3 sin /5+) , (31) 

Y cosh y — sin /3 

I wVcosh^y-sin^/J 

a; = , \6l) 

coshy 

i-fc--<= ""y^':^+^'""^l'y ■ (33) 

cosh y cosh y+ — sin p sin p+ 

In terms of lepton's transverse momentum with respect to the collision axis its rapidity y+ 
and azimuthal angle /3+, the dilepton spectrum reads 

dN^^i+i- _ dN^^i+i- 1 _ e'j^ dN^^i+e- 
dtd'^P-^±dy-^ dtdQ.j^dej^ £+ 7y dt'd^p'_^_ ' 

where (27), (22) where used. The rate on the r.h.s. of (34) is a function of e'^, uj' etc. that should 
be expressed through p+± , /5+ and y+ using (30)-(33) . 

Integration over the lepton direction can be performed if we recall that in the ii''-frame the 
main contribution stems from the coUinear configuration k' ■ v'_^_ 1. Because ycos^+ < 1 and 
l — k'-v'j^ <^ 1, it follows from (24) that 1 — k-V-^- <^ 1. Therefore, integral over the lepton directions 
in K is still dominated by the collinear configuration 6^ ^ 9 and 0+ « (j). Using (22), (19) we get 

= 7ve+(l — cos^ 9) = £_|_sin6' in agreement with (20). To the same approximation dt' de',^ is 
boost-invariant and hence we get in if a formula similar to (12) 



dtd£+ up' 



Ai(Ode+ - Ai'(a;)| , (35) 
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where now 

x = { ^ . ^ ) , x=^sin^. (36) 
V eBs+£- sin 9 J mj 

To calculate the dilcpton spectrum produced by quarks, (34), (35) must be integrated with 
the equivalent photon flux n{uj)duj given by (2). It is helpful to note that the equivalent photon 
spectrum n{u)duj (2) is boost-invariant in our approximation. The differential rate per unit £+, 
y+ and /3+ is rather bulky and we will not write it down explicitly. The rate of lepton production 
with energy Ej^ reads 



TT J,^ a;3 a;[l + (x£/a;)V3] \x ) 

(37) 



dtde^ „ 



where, the invariant parameter x defined in (3) is 



z„eBe smO 



X = ^^^^^^ (38) 
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and ujm is a solution of (4). Integration over £+ as in (15) yields the total lepton rate 

dNg^en- ^ 2a^eBzl r-^du e 2{x^/^ + 1/ Ai' jx) 

dt nmi J^m, xcu u[l + (xeM'/'] J(a/.)V^ xiV4(^3/2 _ 4/^)1/2 ^ • ^ ) 

4. INITIAL QUARK DISTRIBUTION 

To calculate the magnetic contribution to the dilepton spectrum produced in a collision of two 
heavy ions with atomic weights Ai and A2 we need to convolute the initial quark distribution with 
the dilepton spectrum in (37) or (39). For the differential rate per lepton's phase space dr+ we 
have 

Sum runs over all light quarks and antiquarks. In this section we focus on the initial quark 
distribution. 

First, consider dcpendance of the initial quark spectrum on transverse momentum The 
key parameters here arc two saturation momenta Qs\ and Qsi of the two colliding nuclei. They 
depend only on quark's rapidity y and total collision energy. At small quark's transverse momentum 
p_L <C min{Qsi, (5s2}, the spectrum increases as '^^^y^p]^^^ P±) while at p_L S> maxjQsi, (5s2} 
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it falls off as ^^^y^p]^^^ oa l/p\- The maximum is at p± ^ mm{Qsi,Qs2}- It is a reasonable 
approximation to write 



dycPp 

Using this equation in (40) we obtain 

1 



~ 5(p± - mm{Qsi, Qs2}) ] 



2ttp± dy 



(41) 



dtdr. 



2tt ./o J-y/2 dy dtdr 



= E^/ dp I dy 

g 

1 /"^"^ f 

= E - / / 



2- '•^/2 diVA^^,_,gx ciAr^^^+^- 



■ + 



pj_=mm{Qsi,Qs2} 



P±=Qs2 



(42) 
(43) 



where we assumed for definitiveness that Qgi > Qs2 when ?/ > 0. y is the rapidity interval between 
the two ions. 

Initial differential cross section for quark production in heavy-ion collisions in the chiral limit 
reads [8, 14, 17] 



• j d^r'e 



-it-{r-r') ' • 



r ■ r 



daAiA2^qX 4:Nc f 2u [ a2u f ^2 

X ^^Tm { - 'X'- - '■'nSMr - r')) - 1] [SMI - z){r - r'))S2(z{r - r')) - 1] } 



^2 tjrtf'2i 



(r — r'Y" 

X 1 {[5i((l - z)r)S,{zr) - 1] [S^Hl - z)r)S2{zr) - 1]} 



1 

'Zj2 



{ [5i((l - z)r')Si{zr') - l] [52((1 - z)T')S2{zr') - l] } , 



(44) 



where t is the relative transverse momentum of the qq pair, 

Sa{r) = exp I -K'^Qliya, ha) 



(45) 



is the scattering matrix clement, and Qs{ya,ba) is the saturation momentum of nucleus a = 1,2. 
ha is the impact parameter. Integrating (44) over I we obtain the delta function {2TT)'^5{r — r'). 
Subsequent integration over < r < oo and < z < 1 yields 



dcrA-^A2-^qX _ Nc 



J d%^j d%2 



Qil J ^ \ Qs2 / J 



(46) 



dy 247r^ 

It is phenomenologically reasonable approximation to treat the impact parameter dependence 
of Qs as a step function. Denoting by S the overlap area of the two ions, we obtain for the quark 
yield 



dNAiAi-^qX _ NcS 



dy 



247r3 



Qil J V Qs2 / J 



(47) 
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where now Q^'s are taken at ba = 0. Rapidity dependence of the saturation momentum in the 
center-of-mass frame is Qsi,2{y) = Qsoe^^^^^, where A is known phenomenological parameter [18] 
and Qso depends only upon the colhsion energy [19]. With this notation we cast (47) in form 



dNAiA2^qX _ NcSQ 



In (l + e-^^y^ + e'^^ In (l + e^^^^ ] . (48) 



dy 247r3 

Eq. (48) is valid for quark rapidities not too close to the kinematic boundary at y = 
(fragmentation regions of the nuclei). Behavior near the kinematic boundary can be inferred from 
the dependence of the valence quark distribution function qy: when Bjorken's x is close to one 
qv{x) oc (1 — x)^. Using x = (p±/\/s)el^l with Y = ln(s//x^), where s is the center-of-mass energy 
squared and /x Ri 1 GeV we obtain 

dy ~ 247r3 ^^^^ 
where we introduced a shorthand notation 

fiy) = [e^^ln(l + e-2^f) +e-^^ln(l + e2^2')] (^1-^6-^2/1(^/2)^ _ (50) 

5. ELECTRON ENERGY SPECTRUM 

To illustrate the derived results, we numerically compute the magnetic contribution to the 
electron energy spectrum produced in heavy ion collisions at RHIC at the central rapidity y+ = 0. 
Magnetic field strength in this case is such that eB :$> m^. Taking also into account that u ^ me 
we infer from (36) that x S> 1. (This is not the case only at very small angles 9 that are beyond 
the experimental resolution). At large x, the first term in the curly brackets of (35) and (37) is 
negligible compared to the second one. Bearing this in mind and substituting (37) and (48) into 
(43) we derive 



dtde^- Stt^ 
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a^mlY,4 dl3 dyfiy) 
„ -'0 Jo 



X r ^ 1^ / , ,,,,, ( - - >cx'/A Ai'(x) , (51) 

where x, x and % are given by (36), (38) with the following substitutions: 

£ = Qso cosh y e^2//2 ^ sin 6 = y^l - (sin /?/ cosh y)2 . (52) 

Our calculation is valid when £+ ^ e, i.e. in the kinematic region £+ ^ Qgo ~ 1 — 1.5 GeV. For 
electrons produced at j/+ = at RHIC this translates into a condition p+± = e+ < 0.5 GeV. In 
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fact, as we argue in the next section, at larger e+ magnetic contribution to the lepton spectrum is 
negligible. 

Total lepton production rate also simplifies at large x. Taking integral over x in (15) yields [16] 



J (4 



'(4/>.)2/3 xlV4(^3/2 _ 4/^)1/2 

Thus, in place of (39) we have 



2(a;3/2 + l/x)Ai'(x) , ^ o„ 2/3 



0.38- 



J 2m. xV^ u;[l + (xe/u^m ' ^^^^ 



dt ' Trme J2me xVs^; a;[l + (xe/u)^/^] ' 

Convoluting (54) with the hard quark spectrum (49) yields the total electron production rate as a 
function of magnetic field B: 



dt 



da; e 
""L, xV3^'"a;[l + (x£MV3]- ^^^^ 

All equations that we derived so far pertain to lepton production in static magnetic field. In 
fact, magnetic field does change with time albeit adiabatically. Its time-dependence at the central 
rapidity y = is approximately given by [20, 21] 

eB{t) = y exp , (56) 



t2 4t , 

where a is QGP electrical conductivity, Z and Ra are nuclear charge and radius correspondingly. 
Derivation of (56) assumes that QGP conductivity is constant. In expanding medium cr is a 
function of time, e.g. in Bjorken scenario [22] o ~ t"^/^. However it has only mild effect on the 
time-dependence of magnetic field given by (56). Taking for Gold nucleus Z = 79, Ra = 6.5 fm 
we obtain that at t = 0.2 fm (which is ~ l/Qs) ^B = 1.3 in agreement with earlier estimates 
[23]. However, time-dependence of (56) is significantly different from the one found in [23] because 
it takes into account the electromagnetic response of QGP. The time-dependence of magnetic field 
in conducting medium is shown by a solid line in Fig. 3. It is seen that even at t = 10 fm 
magnetic field is about two orders of magnitude larger than the Schwinger field for electrons: 
eBc/m^ = (rTT-e/mTr)^ = 1.3 • 10~^. For comparison, we also calculated the dilepton spectrum 
without the medium effect on magnetic field. In this case, time-dependence of magnetic field can 
be modeled by boosted Coulomb field as 

. 2aZRA7 
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2 4 6 8 10 
tifin) 



FIG. 3: Time-dependence of magnetic field at midrapidity after a collision of two Gold nuclei at y/s~^ = 
200 GeV. Solid line: (56) with cr = 5.8 MeV[24], dashed line: (57). 

where 7 = s/snn f^mN is the Lorentz factor. At RHIC heavy-ions are colhded at 200 GeV per 
nucleon, hence 7 = 100. The time-dependence of magnetic field in vacuum is shown by a dashed 
line in Fig. 3. 

To calculate the total time-integrated yield of electrons substitute (56) (in medium) or (57) 
(in vacuum) into (51) and (55) and integrate over time in the interval 0.2 < t(fm) < 10. The 
result of the calculation is displayed in Fig. 4. Integrating over the electron energy e-|- we obtain 




0.0 0.1 0.2 0.3 0.4 0.5 
s+ (GeV) 



FIG. 4: Spectrum of electrons produced in magnetic field by hard quarks during the QGP life-time. It 
coincides with the distribution in p+±_ at midrapidity y+ = 0. Time dependence of magnetic field is given 
by (56) (sohd hne) and by (57) (dashed line). 

that the total number of leptons produced per event A'^e = 0.049 (in medium), and Ne = 0.0033 
(in vacuum). These numbers as well as Fig. 4 represent a magnetic contribution to e^e~ yield 
at RHIC only due to the square of the amplitude of Fig. 1. It does not take into account other 
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contributions represented by diagrams in Fig. 2. 



6. COMPARISON OF MAGNETIC AND CONVENTIONAL PHOTON DECAY 

MECHANISMS 

In the absence of magnetic field, dilepton production by a fast quark is related to the dilepton 
production by a virtual photon as follows 



drl ~ 3^ / 1^ dT+ 



(58) 



provided that the invariant mass M satisfies M ^> m^. If u; ^ M then the photon spectrum is 
M-independent and can be approximated by that of real photons. In that case we obtain 



dN. 



dV: 



2a oj dN^i 

7^ 111 — ~JF~ ~ "I'^i+t 
37r mi ai + ' 



dN^ 



(59) 



where P, 



describes the probability to produce dilepton via photon decay. This quantity 



should be compared with N^^g+£-, which describes such a probability in magnetic field. If one 
tries to reconstruct the photon spectrum using (58), then the result will be incorrect as it misses 
an important i?-dependent contribution.* To demonstrate how different are these contributions, 
we plotted their ratio in Fig. 5 for azimuthal angle (3 = 0, i.e. perpendicular to B. (We neglect 
small variations with /3). Magnetic contribution dominates at low frequencies. At midrapidity 

100 f 



a. 



10 
1 
0.1 




0.01 



0.0 0.2 0.4 0.6 0.8 1.0 
OJ {GeV) 

FIG. 5: Ratio of dilepton production yields: via photon splitting in magnetic field N^^g+i- and via decay of 
virtual photon without magnetic field Pq^g+i-. Solid and dashed lines correspond to the time-dependence 
of magnetic field according to (56) and (57). Photon's rapidity y = 0. 



OJ 



k± and we can infer from Fig. 5 that at k± = oj < 0.2 GeV the contribution to the dilepton 



* However, if the photon decay vertex is measured, it impiies that photon decay happened long after the disappear- 
ance of QGP and the magnetic field it supported. Such process can certainly be described by (58). 
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yield induced by magnetic field dominates over the conventional virtual photon splitting. Even at 
k± ~ 0.6 GeV it contributes as much as 20% to the yield. Dilepton production in magnetic field 
certainly contributes to resolution of the puzzle of enhancement of dilepton production reported 
in [6]. However, as has been already mentioned, the exact value of magnetic field contribution is 
hard to pin down without accurate knowledge of the magnetic field dynamics. 

7. SUMMARY 

It is hardly surprising that magnetic field generated in heavy ion collisions has a profound impact 
on dilepton production. Indeed, magnetic field strength by far exceeds the critical Schwinger's 
value during the entire QGP lifetime. In the present work we calculated a contribution to dilepton 
spectrum in the region T < e+ < Qg due to magnetic field as a convolution of three factors: 
the initial hard (i.e. non-thermal) quark distribution, equivalent photon fiux and photon decay 
rate. The last two factors exhibit very strong dependence on magnetic field. Because momentum 
is conserved only in the direction of magnetic field, the notion of invariant mass applies only to 
the part of the spectrum independent of magnetic field. Consequently, we plotted the electron 
spectra in Fig. 3 as a function of electron's energy. We also derived formulas for fully differential 
distribution of leptons. We argued that the magnetic contribution is important at electron energies 
below ~ 0.5 GeV at midrapidity at RHIC. In fact it becomes the dominant source of dileptons at 
lower electron energies, see Fig. 5. 

In the region < T contribution of soft quarks, i.e. those quarks that are part of the QGP 
cannot be neglected and must be added to the magnetic contribution of hard quarks computed 
in this paper. Its calculation however is much more complicated because the equivalent photon 
approximation is no longer applicable. We are planning to discuss this contribution elsewhere. 
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